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The theory of mesoscopic fluctuations is applied to inhomogeneous solids consisting of chaotically 
' distributed regions with different crystalline structure. This approach makes it possible to describe 

statistical properties of such mixture by constructing a renormalized Hamiltonian. The relative vol- 
umes occupied by each of the coexisting structures define the corresponding geometric probabilities. 
In the case of a frozen heterophase system these probabilities should be given a priori. And in the 
case of a thermal heterophase mixture the structural probabilities are to be defined self-consistently 
by minimizing a thermodynamical potential. This permits to find the temperature behavior of the 
probabilities which is especially important near the points of structural phase transitions. The pre- 
sense of these structural fluctuations yields a softening of a crystal and a decrease of the effective 
O ' Debye temperature. These effects can be directly seen by nuclear gamma resonance since the oc- 

currence of structural fluctuations is accompanied by a noticeable sagging of the Mossbauer factor 
at the point of structural phase transition. The structural fluctuations also lead to the attenuation 
of sound and increase of isothermic compressibility. 
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I. INTRODUCTION 



There are many examples of matter consisting of regions, chaotically distributed in space, with different structural 
properties. For instance, such are some polymorphic materials. Another example is a crystal subject to strong 
mechanical stress after which the cracks and branches of dislocation are formed in it. These defects have a tendency 
to group inside compact regions. The latter, from the point of view of statistical physics, can be treated as nuclei of the 
amorphised phase inside a crystalline matrix [1]. A similar picture develops in crystals under the action of irradiation 
by fast neutrons when the pores and regions of disorder arise. Under strong irradiation cracks also appear. These 
defects form groups and clusters randomly distributed in space, e.g. as is shown in Figs. 1 and 2. For a statistical 
description of an irradiated crystal the defected regions can be treated as embryos of disordered, usually rarefied, 
£NJ 1 phase inside an ordered, more dense, crystalline structure [2,3]. The relative volume occupied by the disordered phase 
can be measured, with a good accuracy, by investigating the nuclear gamma resonance spectra and the behavior of 
i the Mossbauer factor [4,5]. 

In the considered examples the germs of a disordered structure are randomly distributed in space inside an ordered 
structure. This is why these germs can be called the spatial structural fluctuations. With respect to time, they are 
frozen, which means that their average lifetime, Tf , is much longer than the characteristic time of an experiment, or 
the observation time, r h s , that is: t/ 3> T bs . In the opposite case, when Tf <C r (, s , we have thermal structural 
fluctuations. The example of the latter are even more numerous than those of the frozen structural fluctuations. 

Water, the most widespread matter on Earth, gives us a cogent example of a system consisting of at least two coex- 
isting structures, which is supported by numerous experiments studying its thermodynamic and dielectric properties 
and analysing the Raman spectroscopy data [6-8]. A very convincing is the molecular - dynamic investigation [9] 
showing that fluctuations in water can be decomposed into two components: a fast component ( 1CP 14 — 10~ 13 s ) 
associated with the libration motion in one of the water inherent structures and a slow component ( 10 -12 — 10~ 11 s ) 
associated with the water structure changes. These thermal structural fluctuations in water are related to large local 
energy fluctuations involving about 10-100 molecules. 

Thermal structural fluctuations often appear in the vicinity of structural phase transitions. In principle, these 
fluctuations are possible near both types of structural transitions, near displacive as well as near order-disorder type. 
The distinction between the two types of structural phase transitions can be traced, for example, by their different 
isotope effects [10], although in the presense of strong fluctuations this distinction becomes less pronounced [11 - 
13]. The existence of structural fluctuations around the point of a phase transition reveals itself in the so-called 
pretransitional phenomena [14,15]. 

The freezing - melting phase transition is also a kind of structural transitions. Below the melting point, the 
fluctuational lattice dissociation can arise breaking the local crystalline symmetry, while above the freezing point the 
structural fluctuations are represented by quasicrystalline clusters inside a fluid phase with no crystalline structure. 
Only taking into account the existence of these structural fluctuations it is possible to develop a realistic theory of 
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melting and crystallization [16,17]. Such fluctuations should be especially important during the melting of small 
clusters [18] and solid-fluid phase transition of quantum systems [19]. 

Quantum crystals exhibit near phase transition points the well observed fluctuational coexistence of competing 
structures. This takes place, for example, in solid 4 He along the line of the structural transition between h.c.p. 
and b.c.c. phases and also in the solid mixture of 4 He with 3 He along the stratification line [20]. These local 
fluctuational effects involve usually about 100 particles. 

Fluctuations of one phase inside another can be called [21] the heterophase fluctuations. The structural fluctuations 
are just a sort of the latter. An extensive description of a great number of systems exhibiting these fluctuations has 
been done in the recent review [22]. Therefore, we now shall limit ourselves by the examples considered above, 
mentioning only that the structural fluctuations can play a crucial role in high - temperature superconductors [23] . 

The heterophase fluctuations are mesoscopic since the characteristic size of such a fluctuation, If , satisfies the 
inequality 

a < If < L, 

in which a is the average distance between particles of the system and L is the linear size of the latter. This 
drastically distinguishes the mesoscopic structural fluctuations from the microscopic fluctuations of particles inside 
the same structure. Thus, the oscillations of particles in a crystal in the vicinity of the corresponding lattice sites 
are the most known microscopic fluctuations, the root - mean - square deviation of a particle from its lattice position 
being r <C a . 

Another characteristic feature of the heterophase fluctuations is that the typical energy of each of them is much 
larger than the average single - particle energy [24]. This seems to be quite understandable since a heterophase 
fluctuation, being mesoscopic by its nature, involves many particles whose number is Nf ^ 1 . 

The system with structural fluctuations are difficult to describe. This is because such systems are nonuniform 
with the nonuniformity occurring on mesoscopic scale [22]. Each structural nucleus can have a complicated ramified 
form and can exhibit nontrivial motion [24] . Some thermodynamic features of such systems can be seized by simple 
phenomcnological models (see discussion in Ref.[22] and also [25,26]). However, in phcnomcnological treatment one 
usually has to invoke one or several fitting functions designed to satisfy particular experiments. Of course, this is 
not satisfactory from the point of view of statistical mechanics. These and other difficulties have been discussed and 
illustrated by a mechanical model in Rcf. [27]. 

Recently, a consistent statistical approach has been developed for treating the systems with mesoscopic heterophase 
fluctuations [22]. In the present paper we apply this approach to the systems with structural fluctuations. In Sec. II 
the main steps of constructing a renormalized Hamiltonian for such systems are explained, the general equations 
for the structural probabilities and stability conditions are analysed. In Sec. Ill we concretize the problem for the 
coexisting crystalline structures. The behavior of structure probabilities in the vicinity of a temperature of a struc- 
tural phase transition is studied in Sec. IV, using the Dcbye approximation. A very effective tool for examining the 
properties of structural fluctuations is the Mossbauer effect. Therefore, in Sec.V we scrutinize the peculiar behavior 
of the Mossbauer factor near the temperature of structural phase transition and show how the presense of structural 
fluctuations yields characteristic saggings of the Mossbauer factor. These fluctuations also lead to the attenuation 
of the sound velocity and to the enhancement of the isothermal compressibility, as is discussed in Sec. VI. Resume is 
given in Sec. VII. Everywhere below the system of units is used where Ti = ks = 1 ■ 

II. RENORMALIZED HAMILTONIAN 

Consider a system in which two structures can coexist. We enumerate the structures by the index v = 1, 2 . In the 
case when both structures are crystalline, each of them is characterized by a set 

A„ = {a iu \i = 1,2,... N} (i/ = 1, 2) 

of lattice vectors a iv . As is discused above, the distribution of structures in the space is random. Therefore, it 
is necessary to define the procedure of averaging over structure configurations. Each configuration can be given by 
specifying which regions Y v [y = 1, 2) of the total system volume V are occupied by this or that structure, so that 

V = U„V„ = VxUV 2 . 

Such a specification can be defined [22] by fixing a set 

£ = " = i,2; rev} (l) 
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of functions 



that are called the manifold characteristic functions, or the manifold indicator functions, or simply, the manifold 
indicators. In this way, a distribution of structures in space, that is a structure configuration, is uniquely defined by 
a covering {V^l v = 1,2} of V , or, cquivalently, by the indicator set (1) of manifold indicators (2). The many 
of all possible collections of £ form the topological space T(£) , In order to define the averaging over structure 
configurations, we need to introduce a functional measure on T(£) . This procedure with all mathematical details 
has been thoroughly described in a series of papers [17,28 - 30] and expounded in a recent review [22]. Therefore, 
there is no need to repeat it here. But for the logical self - consistency of the present exposition we will delineate the 
main steps of this procedure. 

Each fixed structure configuration depicts a nonuniform system which can be characterized by the quasiequilibrium 
Gibbs ensemble with a distribution proportional to e~ r (^ , where is a quasi - Hamiltonian defined for a given 

set (1). Specifying a functional measure J D£ on the topological space of structure configurations, we can write the 
partition function as 



-Trf. 



e- r <«Of, (3) 

where Tr means the trace over all quantum - mechanical degrees of freedom, or over the phase space in the classical 
case. The renormalized Hamiltonian H is defined by the relation 

J e- T ®D£ = e-»' T , 

which yields 

H= -Tin / e- T(i) D£,. (4) 



Then the partition function (3) becomes 

Z = Tre- H ' T . 

For each particular structure a space TL U of microscopic states is to be defined, consisting of wave functions 
enjoying the property of the corresponding structural symmetry. The latter, in the case of a crystalline structure, is 
the symmetry of a given crystalline lattice. The renormalized Hamiltonian (4) acts on the fiber space 

y = ® v H v =H l ®H 2 - (5) 

From definition (4) it is clear that the renormalized Hamiltonian should depend on the probability weights of coexisting 
structures. These structure probabilities w v have the meaning of the geometric probabilities, that is, each w v 
determines the ratio of the effective volume occupied by the phase v to the total volume of the system. According 
to this definition, the standard probability properties 

< w v < 1, w v = Wi + U>2 = 1 (6) 

v 

hold. In the case of a frozen structure distribution the structure probabilities w v are to be given a priori. For thermal 
structural fluctuations the values of w v are to be found from the minimization of a thermodynamic potential 

f=-L lnTre -H/T (7) 

under condition (6). Taking into account the normalization condition for w v explicitly, we may write 

Wi = W, W2 = 1 — w. (8) 
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Thence, the extremum of (7) with respect to w is given by the equation 



^ = !<^>=o, 

dw N dw 

where < . . . > implies the statistical average which for an operator A is written as 

1 



(9) 



< A >= TrpA, 



-H/T 



Remind that according to (4) the renormalized Hamiltonian H=H {{w v }) depends on structure probabilities. In 
addition to Eq.(9) defining w v , we have the inequality 



dw 2 



> (0 < w < 1) 



(10) 



showing that the found w u provide us with the minimum of (7). Eq.(10) is the stability condition with respect to 
variations of structure probabilities. From (7) and (10) we have 



1 

N 



d 2 H 1_ [dH 
< dw 2 > T < \ dw 1 ' 



} > 0. 



(11) 



As far as the second term in (11) is always non - negative, the necessary stability condition is 

d 2 u 



< 



dw 2 



> > > 0. 



(12) 



To proceed further, we have to concretize the situation. Denote by Ri v the position of a particle i in a structure 
v and by P iv the momentum of this particle in the same structure; here i = 1,2,...,N and v = 1,2 . Let 
<&{Rijv) be the potential of interaction between the particles i and j for which 



Rijv — Riv Riv 



(13) 



Realizing the procedure described above, after the averaging over structure configurations we obtain [22] the renor- 
malized Hamiltonian 



H= ® V H V = H 1 @H 2 , 



which is a direct sum of terms 



JV 



iV 



2m 2 

2=1 Ijtj 



(14) 



(15) 



where m is a particle mass. 

With the renormalized Hamiltonian given by (14) and (15), Eq.(9) yields 



w = 



2$ 2 + K 2 -K 1 
2($i + $ 2 ) ' 



where 



1=1 



(16) 



(17) 



is the mean kinetic energy per particle in the structure v and 
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is the mean potential energy per particle in the same structure. The necessary condition (12) gives 

$i+$ 2 >0. (19) 

Additional stability conditions follows from the requirement that < w„ < 1 , which leads to the inequalities 

- 2$i <K 1 -K 2 < 2<J> 2 . (20) 

Condition (20) has to be true always since it secures the correct definition (6) of the structure probabilities. If (20) does 
not hold, this means that the considered structures cannot coexist, even in a metastable mixed state, and the system 
is to be in a pure state corresponding to the structure that provides the minimum of the thermodynamic potential 
(7). While when conditions (11), (12) or (19) are not valid but (20) holds, the system with structural fluctuations can 
exist in a metastable state. 

The described approach has been applied to different systems (see [22]). It has made it possible to construct a 
consistent theory of melting and crystalization [16,17] by considering the coexisting solid and liquid structures. It 
is possible to note that when applying this approach to one - and two - dimensional crystals we immediately come 
to the conclusion that in these crystals the infinite long - range order is impossible. Really, suppose that the index 
v = 1 corresponds to a crystalline periodic structure, and v = 2 , to a uniform liquid structure. As is known after 
Peierls [31], the phonon potential energy of one - and two - dimensional crystals diverges. In our case this means that 
<f>i — > oo . Then from (16) it follows that w = wi — > , and, respectively, u>2 — > 1 • This shows that one - and two - 
dimensional periodic structures have zero probability to exist. 

The stability condition (19) has a simple physical interpretation permitting to understand which structures can, in 
principle, coexist. All structures can be roughly divided into two sorts, rigid and soft, according to the sign of . 
The structure can be called rigid if $> v < since this implies that the particles are strongly bound. And we can 
say that the structure is soft if <!>„ > since the particles forming it are weakly bound or unbound. The stability 
condition (19) tells us that two structures can thermally coexist only when at least one of them is soft. Two rigid 
structures cannot thermally coexist in a stable system. For example, liquid has a soft structure. Therefore, a periodic 
crystalline structure and uniform liquid structure can, in principle, coexist near the solid - fluid phase transition point 
[16,17]. Several examples of possible thermal coexistence of different phases have been considered for lattice - gas an 
spin models (see review [22]). Thus, for a lattice - gas model containing two phases with different densities it has been 
shown that such a system is unstable when the effective interaction between particles is attractive [2] but the system 
becomes stable if the effective interaction is repulsive [3]. Of course, lattice - gas models can give only a rough parody 
on real systems with more complicated structures, although these models often describe well chemical processes [32]. 
Low dimensionality (less than three) also makes it more difficult the appearance of thermal structural fluctuations. 
For instance, the latter, as has been rigorously shown [29,33 - 35], do not arise in the two - dimensional Ising - type 
models, though can exist in three - dimensional ones. Nevertheless, a frozen metastable coexistence of two different 
structures seems to be always possible. The condition for the appearance of thermal structural fluctuations changes 
to be more favorable in the vicinity of a structural phase transition where at least one of, or even both, coexisting 
structures become unstable and soft. 



III. CRYSTALLINE STRUCTURES 

Consider the case when two coexisting structures are both crystalline being characterized by different lattice vectors. 

Lattice sites of each structure are given by a vector set { a iv \ i = 1 , 2, . . . , N} . As usual [36] , we can expand $ ( Rij v ) 
in Hamiltonian (15) in powers of the deviations 

from the lattice sites, limiting ourselves by the second order of this expansion. We are interested here in qualitivc 
understanding of the behavior of a heterostructural system, therefore we will not discuss such questions that do 
not change principally this behavior although can be important in quantitative calculations for particular substances. 
These questions include the account of interparticle correlations [17,37 - 40] of anharmonicities [41,42], and of vacancies 

[43]. Instead, we can think of <&(Rij U ) as of an effective potential adjust to take into account these effects, at least 
partially. 
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After expanding Hamiltonian (15) in powers of deviations (21), we get 

H v = E S J + Hf, 

where the first term 

E? = Nw 2 v U v 

is a static potential energy with 

1 N 

Uv =2nH $ (°^) (a l3U =a ll/ - a jv ), 



and the second term 

N 



N 2 2 N 

H ph -w V ^ + ^ V D af3 u a u P 



is the phonon Hamiltonian with the dynamical matrix 



ua iv ua jv 



The mean potential energy (18) becomes 

N .", 

2.Y 



1 N 3 



The Hamiltonian (25) can be diagonalizcd in a standard way by introducing the transformations 



k,s 

• { 9,N J 

k.s 

using the orthogonality condition 

i N 

/ ^ c bi/ c fi:'sV c u kk' u ss' 
i—1 

and defining the frequencies Wfc S! , by the eigenvalue problem 

AT 3 ^ 

^ ^ n a/3tk auB 2 « 



ksv ksv '^hsv* 



771 3 = 1 13=1 

whose eigenfunctions are the polarization vectors e ksu . As a result, we have the phonon Hamiltonian 



Although the way of obtaining (30) is standard, following it we have to be very cautious paying much attention to 
the nontrivial dependence of the Hamiltonian on the structure probability factors w v . Because of this, the phonon 
frequencies defined by (29) become dependent on w v as well. The structure probability w v enters also in all main 
averages such as the phonon distribution 



the momentum squared 



and the correlation function 



^ U+ U ^ \ ( w u^ksu\ -.1 1 

n ksv =< KsMsv >= |exp ^ — - — j-lj , (31) 



iv 2N ^ ^sucoth 2T , (32) 

k.s 



k,s 



In its turn, Eq.(16) for the structure probability involves the mean kinetic energies (17), 

Kv = 47V ^ Uksv 2T ' ^ ' 

k.s 

and the mean potential energies (18), 

$, y = U v + D v , 

^ = ^EEE^<»«.^- (») 

j — 1 k.s a, [3 

The internal energy of the heterostructural system is 

E =<H>= E x +E 2 , 
E v =< H v >= El 1 + Ef, (36) 
where the static energy is given by (23) and the phonon energy is 

Ef ^< H» £ coth ( 37 ) 

k.s 

The latter, using the relations 

K v = w u D Ul w v $ v = w v U v + K v , (38) 

can be written as 

Ef = N{w u K„ + w 2 v D v ) = 2Nw v K v . (39) 

Thus, we see that the internal energy (36) depends on temperature directly and also through the structure probabilities 
w v . Consequently, the specific heat of a heterophase system, 

_ dE ( dE\ , f dE 



8T \dT) w \dw 

contains an additional term, as compared to the specific heat of a pure single - structure system. This excessive term 
makes it possible to explain the so - called specific - heat anomalies observed in heterophase systems [24,44]. 
The free energy (7) takes the form 

l = h + h. 

k.s 

which demonstrates the nonlinear dependence on the structure probabilities w v . 
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IV. STRUCTURE PROBABILITIES 



The equation defining the structure probabilities can be written either as Eq.(16) with substituted there mean 
kinetic energies (34) and potential energies (35) or can be obtained by the direct minimization of (40) with respect to 
w = wi and taking into account (6). Both ways, as can be checked, yield the same answer. To analyse this equation, 
we have first of all to remember that phonon frequencies, given by the eigenvalue problem (29), depend on structure 
probabilities. To make this dependence explicit, we introduce the notation 

ujksv = wl /2 e k su, (41) 



N 3 ^ 

~~ m t—' ^ v ksv ksv 



J=l a,(3 

in which £ksv 

does not contain w v . Emphasizing the dependence on w v explicitly, we have for the kinetic energy 

(34) 

1/2 3/2 
K v = -fij- 2^ s kav coth — . (42) 

k,s 

For the internal energy (36) we get 

E = E l +E 2 , E v =Nw 2 v U v + 2Nw v K v . 

The free energy (40) becomes 

/ = /i + / 2 , 

/, = v£U v + I £ In ^2 sinh . (43) 

Minimizing (43) with respect to w , with the use of notation (8) and relations 

di = dj__di_ df_ 

dw dwi dw 2 ' dw v 

we obtain 



2w v XJ v + 3K V , 



2U 2 + 3(K 2 - Kl ) 

2(U 1 + U 2 ) [ ' 

From the inequalities < w < 1 , assuming that 

Ui + U 2 < (45) 

we have a necessary condition 

U 2 < I (K, - K 2 ) < -Ui (46) 

for w to be considered as a probability. 

To further simplify the analyses, let us resort to the Debye approximation. For this, we pass to the thermodynamic 
limit by using the change 



N 51 n 



' dk 



where 

D = { k I k = I k I < k D } 
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is the Debye sphere, and the Debye momentum ku being defined by the normalization ^ fe g 1 = 1 giving 

k D = (6. 2 p) 1/3 , P^y- 

The phonon spectrum in the Debye approximation acquires the linear form, for which one has to make the substitution 



cl=]imlY( £ -^Y. 
v fc->o 3 ^ V k J 



The quantity c v plays the role of the sound velocity in a pure v - structure. From (41), using the orthogonality 
property 

3 

e ksv e ksv = 

one gets 



N 3 ^2 



In this way, for the kinetic energy (42) we have 

c v = L_ ( 

Jo 



36T 4 /- e -/ 2T 

K » = ™ — / x 3 cothxdx (48) 



and for the free energy /„ we get 



72r 4 r~ 

]v = w z v U v + ——\ xln(2sinhx)dx, (49) 

®v Jo 

where we have introduced the notation 

e„ = w s J 2 T„ D {T vD = c v k D ). (50) 

Here, T^d is the Debye temperature of a pure v - structure, while 9„ can be called an effective Debye temperature 
of a structure inside a mixed hcterostructural system. 

Formulas (48) - (50) show that the low and high temperature limits for a heterostructural system are to be defined 
not with respect to T u d but with respect to the effective temperature Q v given by (50). The latter is renormalized 
by the factor w^ 2 itself depending on temperature. To analyse the behavior of the structural probability (44) we 
need to write accurately the corresponding temperature limits for the kinetic energy (48). 

In the case, when 

r\ 3 / 2 
T<< 2^~ T -< 



4 r B v /2T 



we can use the integrals 



roo 2n-l (2Tr) 2n 



in which B n arc the Bernoulli numbers, 



Bo = 1, Bt = ~, B 2 = 1 B 3 = 0, B 4 = ~, . . . 
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In particular, 



15' 



xdx 



This yields for the kinetic energy (48) 



In the opposite limit, when 



K„ * ±T vDV V* + 3n * T \ /2 . (51) 
!6 lOT^u;, 1 ; 72 



3/2 



2tt 2tt 
using the expansion 

I X 

cothx ~ — I — - — — (x < 7r), 
x 3 45 

we find 

To make the following expressions less cumbersome, it is convenient to introduce the dimensionless static energies 

1 ID 1 ID 



Ux + U 2 



= -(ui+u 2 ). (53) 



-I ID 

In the case of (45), the latter value is positive, u > . Also, we shall use the notation 

-L ID J- ID 

For dcfiniteness, we assume that the structure corresponding to v = 1 has a higher Debye temperature, that is, 
T\d > T 2 d ■ Hence, the parameter t from (54) lies in the region < t < 1 . 

Now, let us understand the behavior of the structural probability (44) in the vicinity of the temperature T s of 
a structural phase transition. This can be a phase transition either between two different crystalline structures or 
between a regular crystalline structure and an irregular glassy structure. The latter consideration is possible owing 
to many similarities between crystalline and glassy states [45,46] and because the Debye approximation is applicable 
to both of them. The qualitative behavior of the structural probabilities near the structural transition temperature 
T s is mainly influenced by the relation between T s and the effective temperatures (50). It is possible to distinguish 
three cases that can be conditionally called the low - temperature, mid - temperature and high - temperature cases. 

Begin with the low - temperature situation when the structural transition temperature satisfies the inequality 

3/2 

T s <^— T vD (i/ = 1,2). (55) 

Then, for the kinetic energies of both structures we can use the approximation (51) which is to be substituted into 
(44). To simplify the resulting expression, we notice that as w u < 1 and T w T s , hence the variable t defined in 
(54) can be considered, according to (55), as small parameter, since t < 1/2-k — 0.159 . This yields 

w ~ ao — ayi 4 , (56) 

where ao is a solution of the equation 

27 _ 27 

uia + — V^o = u 2 (l - a ) + — T\f\ - a , 
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and 

1447t 4 (V^-t 3 VT — ^) 



C(4 



5r 3 [64w v / a (l - a ) - 27(vT 



Remind that w = W\ corresponds to a more rigid structure for which T\£> > T^d ■ Eq.(56) shows that the probability 
of the more rigid structure quickly decreases as temperature increases in the vicinity of the structural - transition 
temperature T s . 

Consider now the case when the transition temperature T s is, in some sense, intermediate satisfying the condition 

3/2 3/2 

^-T 2D <T S < ^-T 1D . (57) 

Then, the kinetic energy of the more rigid structure can be approximated by Eq.(51) while that of the more soft is to 
be approximated by (52), which gives 

9 „ ^_ 3tt 4 T 4 



16 v WT% D y/w 



Substituting this into (44), we find 



where 



w ~ 1 - ait, (58) 



72 

ai 



32ui + 27 

Now again the probability of the more rigid structure decreases with increasing temperature, although not so quickly 
as in (56). The solution (58) exists only if a\ > . If a\ < , we have to put w = 1 , which means that there are 
no structural fluctuations. 

Finally, pass to the high - temperature case, when 

3/2 

T s > -7j- — T v u (i/ = 1,2). (59) 
For the kinetic energies of both structures we can use the approximation (52). Then, (44) yields 

+ £ (60) 



for i> 1 and 



In the vicinity of the structural phase transition the structure probabilities are close to each other, w s=s 1/2 . The 
latter equation becomes asymptotically exact if u\ — > U2 • This is analogous to the case (55) for which (56) also gives 
w w 1/2 if the properties of both coexisting structures are similar to each other, that is if r w 1 and Hi w U2 • 

The above analysis shows that the appearance of thermal structural fluctuations near the point of a structural 
phase transition is facilitated when both coexisting structures have close characteristics. 

V. MOSSB ALTER FACTOR 

The occurrence of structural fluctuations around the point of a structural transition can lead to the emergence 
of various anomalies of observable quantities [22], such as the strong enhancement of specific heat and of diffusion 
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coefficient [22,24,47]. A detailed analysis of experimental data [48 - 51] confirms that, probably, the most common 
feature of structural transitions is an anomalous sagging of the Mossbauer factor near the transition point. Such a 
sagging, as has been proved [52,53], cannot be explained by the existence of a soft phonon mode, which can lead solely 
to a fracture of the Mossbauer factor but by no means to a sagging. However, the general softening of a crystal due 
to the arising structural fluctuations can provoke these saggings, as we demonstrate below. 

The Mossbauer factor of a heterostructural system consisting of two thermally coexisting structures is written in 
the form 

} M = \wnpi + w 2 tp2\ 2 = Jm(T, w), (61) 
where w = w(T) is the structure probability given by (44), and 

ip v = e~ w " = tp v (T,w v ), 

2 2 

W v = q -rl = E q mrl (E q = |-), (62) 



here E q is the recoil energy; q , the absolute value of a gamma - quantum momenta; and 

1 3 

r " = q E < « > ( 63 ) 



3 ^ 



is the mean - square oscillation amplitude of a particle in a v - structure. The correlation function (33) in the Debye 
approximation is 

< u iv u t <W J xcothxdx. (64) 

Whence, the mean - square amplitude (63) becomes 

r l = — 7^r / xcothxdx. (65) 



The low - and high - temperature asymptotes for (65) are 



r 2 



"~ m0 v U 292 J {1 ^ 2^ j ' 

:> ' 7 " H; ' N (T » |^). (66) 



to9„ V©^ v 2tt 

Using (66), and remembering formula (50), for the function ip v defined in (62) we have 



(r>|] 



In order to elucidate the influence of structural fluctuations appearing near the temperature T s of a structural 
transition, we will compare the Mossbauer factor (61) of a heterostructural system with the Mossbauer factor 

rp2 r-T 2D /2T 
^2D JO 



f M (T,0) = exp I -l2E q ^- J xcothxdx} (68) 



of a pure high - temperature structure corresponding to v = 2 . The asymptotic values of the reference factor (68) 



are 
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The change of the Mossbauer factor influenced by the presense of structural fluctuations is convenient to characterize 
by the relative deviation 

xf (r \ _ fM(T,w) 

ojm{T, w) = - 1. (70) 

jmK-L , UJ 

Consider first the case when the temperature of a structural transition is such that 

T s <^ (u = l,2). (71) 

For T w T s we have 

MT,w)*\f M (T,0)] T/2 ^, 



<p 2 (T,w)-[fM(T,0)] 1/2VI ^- 
Therefore, for the relative change (70) we get 



5f M (T s ,w) 



fo 



wfo +( l - w ).fo 



(72) 



where 



f Q = f M (T s ,0), w = w(T s ). 
If we assume that w w 1/2 , then (72) transforms to 



'*y-i. 



When the two coexisting structures are drastically different, so that T\d 3> T 2 d , that is r <C 1 , then 



5/ m (t s ,1)^^(1 + / ™ 7 ) 2 -1 (r«l), 



and when they are similar, so that r w 1 , then 

<5/m(t sj c/r i4 -l (r«l). 

To estimate these quantities, let us take the reference Mossbauer factor f from the region 

0.7 < f < 0.9. 

Then, we obtain 

0.033 < Sf M \ T„ I) < 0.128 (r « 1), 



- 0.137 < 8f M ( T„ -)< -0.043 (r « 1). 



(73) 



(74) 



As we see, a sagging at T s can be directed upward as well as downward, depending on the parameter t = TidITid 
characterizing the difference between the structures. When t ss 0.5 , there is no sagging at all. 
Turn now to the case when the temperature of a structural transition is in the interval 



2n <ls< 2tt- 



(75) 
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Hence, at T w T s one gets 



w (r,»).[/ M (r,o)r ! / 8 ^ 



^(T,l-«;)~[/ M (T,0)] 



l/2(l-tu) 2 



For the relative change (70) we find 



8fM{T s ,w) 



f 



o 



w 2 /8t sy /w 



(1 - H/o 



1/2(1-k,) 2 



- 1 



When T s is in the middle of the interval (75), we can use the approximation 



-t ID £ 



Taking w w 1/2 , we reduce (76) to 



^( T 4)"4k( / o 2/2 ^ (1+T)+/ o) 2 



1. 



Considering again two limiting situations of very different and similar structures, we 

£/m(t 5 ,0 ~-L(l + /2) 2 -l (r«l), 



2/ 4/ 



<5/ M (T s ,^)^ — (/ °- 177 + / 2 ) 2 -l (r«l). 



Taking the values of /o from (73), we come to the result 



1 



-0.207 < 5f M ( T s , - ) < -0.090 (r « 1), 



0.271 < 5f M [ T a ,- )< -0.108 (twI). 



Here the sagging are always directed downward. 

If the structural transition occurs at high temperature, when 



then 



Ts> 2^ (^ = 1 ' 2 )' 



^(r, w )c.[/ M (r,o)r 2 /^ 2 , 



^(t.^-i/mCt.o)] 1 /^-™) 3 

at T «T S . The relative change (70) becomes 



For u; w 1/2 Eq. (79) gives 



Sf M (T a , W )*T wfl 2/2w2 + (l- W )fl /2il - w) ' 
Jo 



- 1. 
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from where 



<\/w (tJ-^j ^-L(i + f*) 2 -i (r «i), 



Thence, invoking (73), we obtain 



-0.207 < 5/ M (t s , < -0.090 (t « 1), 

- 0.657 < <5/m (t s , < -0.271 (r w 1). (80) 

All saggings are again directed downward. 

Thus, we see that the appearance of structural fluctuations near the temperature of a structural phase transition can 
yield a noticeable sagging of the Mossbauer factor as a function of temperature. In the majority of cases this sagging 
is directed downwards. The sagging can be easily observed if usil/2, and immediately disappears when any of the 
structure probabilities tends to zero. For example, suppose that the probability of the structure fluctuations of a high 
- temperature more soft phase inside a low - temperature more rigid structure is very small, so that w 2 = x <C 1 . 
Expanding (61) in powers of x , we have 



/m(T,1-z)~/ m (T,1) 



2 + ^ 



2Ti 



D 



In the standard Mossbauer experiments the recoil energy is negligibly small as compared to the Debye temperature, 
E q <C Tin , usually E q /T\B ~ 10~ 7 . Therefore, the change of the Mossbauer factor 

f M {T,l-x)~f M {T,l){l-2x) (x«l) 

is also quite small and disappears as soon as x — > . 

To illustrate that the values of saggings at T s are in agreement with experiment, let us consider the structural 
transition between the low - temperature cubic phase and the high - temperature rhombic phase in the compounds 
Sn\- x Ge x Te . The characteristic temperatures of the latter are T s w and T\r> ~ T 2 d ~ 150K . Mossbauer 

investigation [54] display the existence at T s of a pronounced sagging of the Mossbauer factor, Sf^ p (T s ) w —0.4 . 
The considered characteristic temperatures correspond to inequality (78) and to t m 1 . The measured Mossbauer 
sagging is in agreement with the second line of estimate (80). 

VI. SOUND VELOCITY 

The existence of structural fluctuations near a phase transition point can also lead to a distinct attenuation of the 
velocity of sound 



lim y 

k^0 k ^ 



' < Hf > . (81) 



8n ksv 



To trace out the dependence of the sound velocity (81) on the structure probabilities, let us introduce the notation 

c sv = lim ^ (82) 

k^O k 

for the sound velocity with polarization s inside a pure structure v . According to (30) and (31), the sound velocity 
(81) in a heterostructural substance takes the form 

3/2 . 3/2 , 00 n 

v a = w 1 ' c s i+w 2 ' c s2 . (83) 
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Defining the average, with respect to polarizations, velocities of sound for a heterostructural system, 



1 3 



(84) 



and for a pure v - structure, 



we get 



1 3 

"v = 77 y ' c sv ! 



3/2 . 3/2 
V = W{ Ci + 10 2 ' C2, 



(85) 



where each of c„ is given by (47). 

The relative decrease of the sound velocity, due to arising structural fluctuations, can be characterized by 



In the Debye approximation we have 



Therefore, (86) gives 



8v(w) = 1. 

C2 



C2 T 2 D 



(86) 



Cl 



- ID 



5v(w) = -w 3 / 2 + (1 — w) 3 / 2 — 1. 



(87) 



Assuming that at the temperature T s of a structural transition one has w w 1/2 and r ~ 1 , we obtain 



Sv l^-j w -0.293 (T = T s ). 

This decrease of the sound velocity is completely due to the onset of mesoscopic structural fluctuations and can happen 
even at first - order phase transitions when there are no microscopic critical fluctuations [55] related to second - order 
phase transitions. For example, a similar decrease of the sound velocity occurres at freezing point of water [56]. 
As far as the average sound velocity (84) can be expressed through the derivative 



dP 



mv 



of pressure with respect to density, and the same derivative is involved into the definition of the isothermal compress- 
ibility 



Kt = -[V 



,dP_ 
dV 



= \ P 



opy 1 
'dp) 



we can easily find the influence of structural fluctuations on the latter. Thus, introducing the isothermal compressibility 
of pure structures, 



1 



mpcl ' 



for a heterostructural system we find 



Kt 



3\ V2 



3\ 1/2' 



K 2 



(88) 
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Defining the relative change of compressibility by 



5kt(w) = 1 



(89) 



and using the previous notation, we come to 



5kt{w) 



[ w 3/2 _|_ T (i 




(90) 



If at the point of a structural phase transition we have w 



1/2 and r w 1 , then 



S*r (i) « 1 



(T = T S ). 



This means a strong increase of compressibility. 



VII. CONCLUSION 



The main aim of this paper has been to present a general approach for describing statistical properties of heterostruc- 
tural systems. This approach can be used for treating heterogeneous system with frozen structural fluctuations induced 
e.g. by shock waves or irradiation. Then, the relative volumes occupied by each structure, or the structure proba- 
bilities, are additional thermodynamic variables which can be defined experimentally by different means, for example 
by nuclear gamma resonance [4,5]. Probably, the most promising application of this approach is to considering sub- 
stances with thermal structural fluctuations. This is of great importance for describing systems with structural phase 
transitions. The appearance of structural fluctuations around the phase transition point leads to various pretransi- 
tional, or precursor, phenomena that are often manifested in pronounced anomalies of thermodynamic and dynamic 
characteristics. The liquid - solid phase transition can be regarded as a kind of such structural transitions [16,17]. 
There are other numerous examples of structural transitions accompanied by the occurrence of structural fluctuations 
whose existence can be observed by different experiments. Some of the examples have been discussed in this paper. 

A number of other examples has been reviewed recently [22] . As is well known, structural fluctuations play a decisive 
role in high - temperature superconductors [23], in solids with martensitic transformations [57,58], and in crystals 
with perovskite structure whose structural transitions are accompanied by experimentally well observed pretransitional 
structural fluctuations [59,60]. A more detailed theoretical consideration of these particular substances is supposed 
to be done in separate publications, basing on the approach developed in the present paper. 
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Figure Captions 

Fig.l. The typical distribution of defect clusters inside a metal irradiated by fast neutrons. 

Fig. 2. Pores and cracks in a metal irradiated by fast neutrons: darker regions correspond to the amorphised phase 
with disordered structure. 
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